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Abstract These lectures give an introduction to the theory of holographic super- 
conductors. These are superconductors that have a dual gravitational description 
using gauge/gravity duality. After introducing a suitable gravitational theory, we 
discuss its properties in various regiemes: the probe limit, the effects of backre- 
action, the zero temperature limit, and the addition of magnetic fields. Using the 
gauge/gravity dictionary, these properties reproduce many of the standard features 
of superconductors. Some familiarity with gauge/gravity duality is assumed. A list 
of open problems is included at the end. 



1 Introduction 



The name "holographic superconductor" suggests that one can look at a two (spa- 
tial) dimensional superconductor and see a three dimensional image. We will see 
that there is a class of superconductors for which this is true, but the image one 
"sees" is quite striking. It involves a charged black hole with nontrivial "hair". This 
remarkable connection between condensed matter and gravitational physics was dis- 
covered just a few years ago. It grew out of the gauge/gravity duality which has 
emerged from string theory ll47ll2Tll6T1l . Although this duality was first formulated 
as a equivalence between a certain gauge theory and a theory of quantum grav- 
ity (and provided new insights into each of these theories), over the past decade 
it has been applied with notable success to other areas of physics as well. Many 
of these new applications are discussed in other lectures in this school. I will fo- 
cus on the application to superconductivity. These lectures will be heavily based on 
ll29l [30l [37 38 1 . For a more general discussion of applying gauge/gravity duality 
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to condensed matter, see the excellent reviews by Hartnoll ll27l . Herzog (32 1 and 
McGreevy||48l. 

I will start with a brief introduction to superconductivity. In section 2, I will 
introduce a simple model for a holographic superconductor. Most of our discussion 
will be devoted to exploring the consequences of this model, beginning in section 
3 with the probe limit - a simplification of the model which preserves most of the 
physics. In section 4 we will discuss the full solution with backreaction. We next 
examine the ground state of the superconductor (section 5), and study its behavior 
when magnetic fields are added (section 6). We conclude with a brief discussion of 
recent developments (section 7) and a list of open problems (section 8). 



1.1 Superconductivity 

It was noticed in the early part of the 20' /! century that the electrical resistivity of 
most metals drops suddenly to zero as the temperature is lowered below a critical 
temperature T c . These materials were called superconductor^] A second indepen- 
dent property of these materials was the Meissner effect: A magnetic field is ex- 
pelled when T < T c . This is perfect diamagnetism and does not follow from the 
perfect conductivity (which alone would imply that a pre-existing magnetic field is 
trapped inside the sample). A phenomenological description of both of these proper- 
ties was first given by the London brothers in 1935 with the simple equation oc A, 
B4l . Taking a time derivative yields Ej dJ//dt, showing that electric fields ac- 
celerate superconducting electrons rather than keeping their velocity constant as in 
Ohm's law with finite conductivity. Taking the curl of both sides and combining 
with Maxwell's equations yields V 2 Bj °< Bj showing the decay of magnetic fields 
inside a superconductor. 

In 1950, Landau and Ginzburg described superconductivity in terms of a second 
order phase transition whose order parameter is a complex scalar field <p [ 17|. The 
density of superconducting electrons is given by n s = |<p(x)| 2 . The contribution of 
(p to the free energy is assumed to take the form 

F = a(T-T c )\<p\ 2 + ^\<p\ 4 + --- (1) 

where a and j3 are positive constants and the dots denote gradient terms and higher 
powers of (p. Clearly for T > T c the minimum of the free energy is at <p = while 
for T < T c the minimum is at a nonzero value of (p. This is just like the Higgs 
mechanism in particle physics, and is associated with breaking a 1/(1) symmetry. 
The London equation follows from this spontaneous symmetry breaking |60|. 

A more complete theory of superconductivity was given by Bardeen, Cooper 
and Schrieffer in 1957 and is known as BCS theory 0. They showed that interac- 
tions with phonons can cause pairs of elections with opposite spin to bind and form 
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a charged boson called a Cooper pair. Below a critical temperature T c , there is a 
second order phase transition and these bosons condense. The DC conductivity be- 
comes infinite producing a superconductor. The pairs are are not bound very tightly 
and typically have a size which is much larger than the lattice spacing. In the super- 
conducting ground state, there is an energy gap A for charged excitations. This gap 
is typically related to the critical temperature by A ss 1 HT C . The charged excitations 
are "dressed electrons" called quasiparticles. The gap in the spectrum results in a 
gap in the (frequency dependent) optical conductivity. If a photon of frequency co 
hits the superconductor, it must produce two quasiparticles. The binding energy of 
the Cooper pair is very small, but the energy of each quasiparticle is A , so the gap 
in the optical conductivity is co g = 2A s» 3.5T C . 

It was once thought that the highest T c for a BCS superconductor was around 
30" K. But in 2001, MgZ?2 was found to be superconducting at 4-0" K and is believed 
to be described by BCS. Some people now speculate that BCS could describe a 
superconductor with T c = 200" Kr\ 

A new class of high T c superconductors were discovered in 1986 17] ■ They are 
cuprates and the superconductivity is along the CuOj. planes. The highest T c known 
today (at atmospheric pressure) is T c — 134" K for a mercury, barium, copper oxide 
compound. If you apply pressure, T c climbs to about 160K. Another class of super- 
conductors were discovered in 2008 based on iron and not copper l40l . The highest 
T c so far is 56" K. These materials are also layered and the superconductivity is again 
associated with the two dimensional planes. They are called iron pnictides since they 
involve other elements like arsenic in the nitrogen group of the periodic table. 

There is evidence that electron pairs still form in these high T c materials, but the 
pairing mechanism is not well understood. Unlike BCS theory, it involves strong 
coupling. Gauge/gravity duality is an new tool to study strongly coupled field the- 
ories. In particular, it allows one to compute dynamical transport properties of 
strongly coupled systems at nonzero temperature. Condensed matter theorists have 
very few tools to do this. We will describe below the first steps toward applying this 
new tool to better understand high T c superconductivity. I must stress at the begin- 
ning that we are still at the early stages of this endeavor. We will construct simple 
gravity models and show that they reproduce basic properties of superconductors. 
But our models are too crude to make detailed comparisons with any real-world 
material. 



2 A gravitational dual 

How do we go about constructing a holographic dual for a superconductor? The 
minimal ingredients are the following. In the superconductor we need a notion of 
temperature. On the gravity side, that role is played by a black hole. Recall that in 
the 1970's Hawking (following work by Bekenstein and others) showed that sta- 
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tionary black holes are thermodynamic objects with a temperature T related to the 
surface gravity K via T — k/2tc. In gauge/gravity duality, the Hawking temperature 
of the black hole is identified with the temperature of the dual field theoruj Since 
gauge/gravity duality traditionally requires that spacetime asymptotically approach 
anti de Sitter (AdS) space at infinity, we will be studying black holes in AdS. Un- 
like asymptotically flat black holes, these black holes have the property that at large 
radius, their temperature increases with their mass, i.e., they have positive specific 
heat, just like familiar nongravitational systems. There are also planar AdS black 
holes, which will be of most interest. These black holes always have positive spe- 
cific heat. 

In the superconductor, we also need a condensate. In the bulk, this is described by 
some field coupled to gravity. A nonzero condensate corresponds to a static nonzero 
field outside a black hole. This is usually called black hole "hair". So to describe 
a superconductor, we need to find a black hole that has hair at low temperatures, 
but no hair at high temperatures. More precisely, we need the usual Schwarzschild 
or Reissner-Nordstrom AdS black hole (which exists for all temperatures) to be un- 
stable to forming hair at low temperature. At first sight, this is not an easy task. 
There are "no-hair theorems" which say that certain matter fields must be trivial 
outside a black hole (see, e.g., |8] [35)). The idea behind these theorems is simply 
that matter outside a black hole wants to fall into the horizon (or radiate out to 
infinity in the asymptotically flat case). However, there is no general "no-hair the- 
orem". Each matter field must be considered separately. The result is a set of black 
hole uniqueness theorems showing that when gravity is coupled to certain matter 
fields, stationary black holes are uniquely characterized by their conserved charges: 
mass, angular momentum and electromagnetic charge. These theorems usually re- 
quire linear matter fields or scalars with certain potentials V((j>). Counterexamples 
to a general no-hair theorem have been known since the early 1990's. (So our task 
is not impossible.) For example, it was shown that static Yang-Mills fields can exist 
outside the horizon |59|. 

String theory has many dilatonic black holes with scalar hair, and one might 
be tempted to try to use one of these to model a superconductor. But if the action 
includes a term like e 2 "^ F^F^ ', this is doomed to failure. In this case, F 2 is a 
source for <p, so all charged black holes have nonzero <$>. This "secondary hair" is 
not what we want, since we want the hair to go away at high temperatures. (Theories 
with more general coupling f(ij))F^ v F^ v are possible candidates provided / does 
not have a linear term in (j). However the example we will study uses a standard 
Maxwell action.) A general argument against AdS black holes developing scalar hair 
was given by Hertog [31 1. He considered a real scalar field with arbitrary potential 
V((j>) (with negative extremum so AdS is a solution), and showed that neutral AdS 
black holes can have scalar hair if and only if AdS itself is unstable. This is clearly 
unacceptable. 

A surprisingly simple solution to this problem was found by Gubser |20l . He 
argued that a charged scalar field around a charged black hole in AdS would have 



3 I will assume that the reader is familiar with the basics of gauge/gravity duality. If not, see other 
contributions in this book. 
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the desired property. Consider 

S = j 'd A x^~g[R+ 1-^ V F^ V -|Vf -iqAV\ 2 -m 2 |f 2 |) . (2) 

This is just general relativity with a negative cosmological constant A = —3/L 2 , 
coupled to a Maxwell field and charged scalar with mass m and charge q. It is easy 
to see why black holes in this theory might be unstable to forming scalar hair: For 
an electrically charged black hole, the effective mass of f is m 2 ff — m 2 + q 2 g"A 2 . 
But the last term is negative, so there is a chance that m 2 ff becomes sufficiently neg- 
ative near the horizon to destabilize the scalar field. Furthermore, as one lowers the 
temperature of a charged black hole, it becomes closer to extremality which means 
that g„ is closer to developing a double zero at the horizon. This means that \g"\ 
becomes larger and the potential instability becomes stronger at low temperature. 

We will see that black holes in this theory indeed develop scalar hair at low 
temperature. One might wonder why such a simple type of hair was not noticed 
earlier. One reason is that this does not work for asymptotically flat black holes. 
The AdS boundary conditions are crucial. One way to understand the difference is 
by the following quantum argument. Let Qj be the initial charge on the black hole. 
If qQi is large enough, even maximally charged black holes with zero Hawking 
temperature create pairs of charged particles. This is simply due to the fact that the 
electric field near the horizon is strong enough to pull pairs of oppositely charged 
particles out of the vacuum via the Schwinger mechanism of ordinary field theory. 
The particle with opposite charge to the black hole falls into the horizon, reducing 
Qi while the particle with the same sign charge as the black hole is repelled away. 
In asymptotically flat spacetime, these particles escape to infinity, so the final result 
is a standard Reissner-Nordstrom black hole with final charge Qf < Qi- In AdS, the 
charged particles cannot escape since the negative cosmological constant acts like a 
confining box, and they settle outside the horizon. This gas of charged particles is 
the quantum description of the hair. This quantum process has an entirely classical 
analog in terms of superradiance of the charged scalar field. 

The four dimensional bulk theory |2]) is dual to a 2+1 dimensional boundary the- 
ory. This is the right context to try to understand the superconductivity associated 
with two dimensional planes in the high T c cuprates or iron pnictides. One can also 
study this theory in five dimensions to describe 3 + 1 dimensional superconductors. 
I should emphasize that at the moment we are not trying to derive the gravitational 
theory from string theory. The idea is to find a simple gravity theory with the prop- 
erties we want, and analyze it using standard entries in the gauge/gravity duality 
dictionary. However, we will see later that this simple model can, in fact, be realized 
as a consistent truncation of string theory. 

Before we proceed, I should comment on the following confusing point. In 
gauge/gravity duality, gauge symmetries in the bulk correspond to global symme- 
tries in the dual field theory. So although the scalar hair breaks a local 1/(1) symme- 
try in the bulk, the dual description consists of a condensate breaking a global t/(l) 
symmetry. Thus strictly speaking, the dual theory is a superfluid rather than a su- 
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perconductor. The superfluid properties of this model have indeed been investigated 
ll34l |6l [3ll. However, one can still view the dual theory as a superconductor in the 
limit that the U(l ) symmetry is "weakly gauged']^] In fact, most of condensed mat- 
ter physics does not include dynamical photons, since their effects are usually small. 
In particular, BCS theory only includes the electrons and phonons. Electromagnetic 
fields are usually introduced as external fields, as we will do here. (Unfortunately, 
our electromagnetic fields will not be dynamical on the boundary.) 



3 Probe Limit 

If one rescales = A^/q and W = W/q, then the matter action in (|2ji has a 1/q 2 in 
front, so the backreaction of the matter fields on the metric is suppressed when q is 
large. The limit g — !• °° with A^ and ¥ fixed is called the probe limit. It simplifies the 
problem but retains most of the interesting physics since the nonlinear interactions 
between the scalar and Maxwell field are retained. In this section, we will explore 
this probe limit. We first discuss the formation of the condensate and then compute 
the conductivity. To simplify the presentation, we will drop the tildes. 



3.1 Condensate 



We start with the planar Schwarzschild anti-de Sitter black hole in four dimensions 

dr 2 

ds 2 = -f{r)dt 2 + j^+r 2 {dx 2 +dy 2 ), (3) 



where 



2 / J3 



r I . r. 



/= zn 1_ #' ■ (4) 



L is the AdS radius, and the Schwarzschild radius rrj determines the Hawking tem- 
perature of the black hole: 

T =4th- (5) 

In the probe approximation this metric is a fixed background in which we solve the 
Maxwell-scalar equations. 

We assume a plane symmetric ansatz, 

W=V(r), A t = <t>(r) (6) 



4 This means that one imagines that the dual action includes terms like | V, — ieAi)(p\ 2 with very 
small charge e. 
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With A r = A x = A y = 0, the Maxwell equations imply that the phase of \jf must be 
constant. Without loss of generality we therefore take \j/ to be real. The Maxwell- 
scalar field equations reduce to the following coupled, nonlinear, ordinary differen- 
tial equations: 

7' 2\ , 2 m 2 



V + ( - + z ) V +72V-— V = ( ? ) 



„ 2 , 2w 2 
f +-f --^</>=0. (8) 

The key term in the first equation is (0 2 // 2 ) y/. This comes in with the opposite sign 
of the mass term and will cause the scalar hair to form at low temperature. 

We first consider the case m 2 = —2/L 2 . It might seem strange to make the scalar 
field tachyonic, but this mass is perfectly allowed in gauge/gravity duality. First note 
that a tachyonic mass in field theory does not describe particles moving faster than 
the speed of light. Instead, it usually describes an instability. The value y = is 
unstable and the field rolls off the potential. However, Breitenlohner and Freedman 
151 showed that AdSd+\ spacetime is stable even with scalar fields with m 2 < 
provided m 2 > mjj F with 

, d 2 

m BF = -^2 (9) 

This is because there is so much volume at large radius that the positive gradient 
energy can compensate for a negative m 2 . Our choice of mass satisfies this bound 
and in fact corresponds to a conformally coupled scalar in AdS$. The standard com- 
pactification of 11D supergravity on S 7 produces many fields with this mass. 

We now consider the boundary conditions. At the horizon, one often argues that 
(j) = A t must vanish in order for g^A^Ay to remain finite. This is correct, but A^ is 
gauge dependent, so it is not obvious that a diverging vector potential is a problem if 
the Maxwell field remains finite. A better argument for setting A, = at the horizon 
is the following The source for Maxwell's equations in the bulk is, of course, 
gauge invariant. But in a gauge in which y/ is real, the current is just y/ 2 A^ . Since the 
current must remain finite at the horizon, we need A^ to remain finite, and hence = 
A, = 0. Even in Einstein-Maxwell theory (without charged sources) A, must vanish 
at a static black hole horizon by the following argument: For describing thermal 
properties of the black hole, one should use the Euclidean solution. The Wilson 
loop of A^ around the Euclidean time circle is finite and gauge invariant. If A, is 
nonzero at the horizon, the Wilson loop is nonzero around a vanishing circle which 
implies that the Maxwell field is singular. 

There is another constraint that must be satisfied in order for the solution to 
be smooth at the horizon. Mutliplying (j7]i by / and evaluating at r = one finds 
f'\j/' — m 2 \j/, so i/A(ro) and V / '( r o) are not independent. As a result, even though 
we start with two second order equations which have a four parameter family of 
solutions, there is only a two parameter subfamily which is regular at the horizon. 



5 I thank Karl Landsteiner for suggesting this. 
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They can be labelled by y/(ro) and (j)'(ro). (Note that <j>'(ro) is essentially the electric 
field at the horizon: ((j)') 2 = F^ V F» V .) 

We now turn to the boundary condition at infinity. Asymptotically: 

w =^— + ^- + .... (10) 

and 

</)=Ai-- + -"- (ID 
r 

Usually, normalizability requires that the leading coefficient in y/ must vanish. How- 
ever, since we have chosen a mass close the the BF bound, even the leading term in 
\j/ is normalizable. In this case, one has a choice: One can consider solutions with 
yW — o or y/ 2 ) = 0. For definiteness, we will mostly consider standard boundary 
conditions, y/ 1 ) = . 

After imposing this asymptotic boundary condition, we have a one parameter 
family of solutions, which can be found numerically 1 29 1 . (See |19] for an accurate 
analytic approximation and [43 1 for an analytic solution of a related model.) We 
won't present any plots of y/(r) and (j)(r) since they look rather boring. It is easy 
to see from (|8]l that (j) is a monotonic function. It starts at zero at the horizon and 
at any local extremum (j)" °< 0. So it cannot have a positive maximum or a negative 
minimum. If it starts increasing away from the horizon, it continues to increase and 
asymptotically approaches the constant jj,. y/(r) does not have to be monotonic. 
There is a discrete infinite family of solutions for y/(r) that satisfy our asymptotic 
boundary conditions. They can be labelled by the number of times y/ vanishes. It is 
believed that only the lowest solution which monotically decreases from y/(ro) to 
zero is stable, although I do not think a stability analysis has been performed yet. 
Despite the fact that the solutions of interest are simple monotonic functions, they 
have important consequences for the dual field theory as we now discuss. 

Since we have not started with a consistent truncation of string theory, we do 
not have a detailed microscopic description of the dual field theory. Nevertheless, 
basic elements of the gauge/gravity duality dictionary allow us to say the following. 
The dual theory is a 2+1 dimensional conformal field theory (CFT) at temperature 
T given by |5j. The local gauge symmetry in the bulk corresponds to a global 1/(1) 
symmetry in the CFT. The asymptotic behavior of the bulk solution determines cer- 



tain properties of the dual field theory. For example, from (111, jU. is the chemical 
potential and p is the charge density. It may seem strange that our superconductor 
has a nonzero charge density. After all, ordinary superconductors are electrically 
neutral. Perhaps the best analogy is to say that our holographic superconductor is 
modeling the electrons, but does not include the (positively charged) atomic lattice. 
Indeed, our model has translational symmetry, and there is no sign of any lattice. 
From a practical standpoint, we need the charge since neutral black holes in the bulk 
are not unstable to forming scalar hair. In addition, without the charge (or chemical 
potential), the dual theory is scale invariant and cannot have a phase transition. 

The dual theory also has an operator charged under the t/(l) which is dual to y/. 
Since we have chosen m 2 close to the BF bound, there are two possible operators 
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depending on how one quantizes y/ in the bullae 11421 , If the modes are defined with 
the standard boundary condition (faster falloff) for y/ in the bulk, the dual operator 
has dimension two. In this case, a nonzero y/ 1 ) corresponds to adding a source for 
this operator in the CFT, and a nonzero y/ 2 ' corresponds to a nonzero expectation 
valuqj 

o 2 = y {2) . (12) 

Since we want the condensate to turn on without being sourced, we have set y/ 1 ) = 
0. There is an alternative quantization of y/ in the bulk in which the roles of y/ 1 ' 
and y/ 2 ' are reversed. y/ is now dual to a dimension one operator. If one wishes to 
study this case, one should impose the boundary condition y/ 2 ' = 0. 

We want to know how the condensate O2 behaves as a function of temperature. 
Before presenting the results we have to discuss an important scaling symmetry of 
our problem. In any conformal field theory on R", one can change the temperature 
by a simple rescaling. In the bulk, this is reflected in the statement that the rescaling 

r — > ar, (t,x,y) — > (t,x,y)/a, ro — > aro (13) 

leaves the form of the black hole ^ invariant with / — > a 2 f. It is easy to check 
that the Maxwell-scalar field equations ( |7|8| l are invariant under this rescaling if 
<p — > a<p (so A = (j)dt is invariant) and y/ is unchanged y/ — » y/. Rather than discuss 
this trivial dependence on temperature which simply reflects the scaling dimension, 
we are interested in the dependence of a dimensionless measure of the condensate 
as a function of a dimensionless measure of the temperature. It is convenient to use 
the chemical potential to fix a scale and consider ^/Ch,/^ as a function of T/fi. 
When one does this, one finds that the condensate is nonzero only when T/jj, is 
small enough. Setting T c equal to the critical temperature when the condensate first 
turns on, we get Fig. 1. 

This curve is qualitatively similar to that obtained in BCS theory, and observed 
in many materials, where the condensate rises quickly as the system is cooled below 
the critical temperature and goes to a constant as T — >• 0. Near T c , there is a square 
root behavior O2 = 100r ( 2 (l — T /T c ) l l 2 . This is the standard behavior predicted by 
Landau-Ginzburg theory. 

A nonzero condensate means that the black hole in the bulk has developed scalar 
hair. One can compute the free energy (euclidean action) of these hairy configu- 
rations and compare with the solution y/ = 0, <p = p(l/Vo— l/r) with describes a 
black hole with the same charge or chemical potential, but no scalar hair. It turns 
out that the free energy is always lower for the hairy configurations and becomes 
equal as T — >• T c ||29ll . The difference of free energies scales like (T c — T) 2 near the 
transition, showing that this is a second order phase transition. Actually, you can 
compare the free energy at fixed charge or fixed chemical potential. These corre- 
spond to two different ensembles. In both cases, the free energy is lower for the 
hairy configuration. 



6 We will work in the classical limit, but the correspondence applies to the full quantum theory. 

7 This normalization of O differs from that of 1 29 1 by a factor of v2. 
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Fig. 1 The condensate as a function of temperature. The critical temperature is proportional to the 
chemical potential. 



We now generalize to other masses. Recall that the asymptotic behavior of a 
scalar field of mass m in A1JS4 is 

yA— yAj^ 

where 



A± = i^3±^9+4(mL)2j. (15) 

For m 2 > m 2 BF + L~ 2 , only the mode with A + falloff is normalizable, and so we 
interperet y/ + = O, where O is the expectation value of the operator dual to the 
scalar field, and \ + is the dimension of the operator. is dual to a source for 
this operator, so we only consider solutions with the standard boundary condition 
= 0. For mg F < m 2 < m\ F +L~ 2 one can also consider the alternative boundary 
condition y/ + = 0, and the dual operator has dimension A_ . 

Fig. 2 shows the results for the condensate for various masses. It is convenient to 
label the curves in terms of the dimension of the condensate A, rather than the mass 
to distinguish the two possible cases when m 2 = —2/L 2 . The masses range from the 
BF bound to m 2 = 0. Similar behavior is found for m 2 > BTIl . 

The qualitative behavior is the same as before. In all cases, there is a critical tem- 
perature T c (proportional to the chemical potential) above which the condensate is 
zero. Near the critical temperature, O °^ (T c . — T) l l 2 . In all but one case, the con- 
densate saturates as T — > 0. The exceptional case is the dimension one curve which 
starts to grow at low temperature. When the condensate becomes very large, the 
backreaction on the bulk metric can no longer be neglected. We will see later that in 
the full solution with backreaction, the condensate approaches a finite limit at zero 
temperature. 
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Fig. 2 Condensates with different dimension, X, as a function of temperature. The condensate 
tends to increase with X. Figure is taken from |37|. 



It turns out that the ratio T c /ji increases as the dimension of the operator de- 
creases. For A > 3/2, one can understand this since a smaller A corresponds to a 
smaller m 2 making is easier for the scalar hair to form. However, this continues for 
A < 3/2 when one must increase m 2 (and use the alternative boundary conditions) 
to decrease the dimension. In fact, as observed in [ 12], T c / ji appears to diverge as A 
approaches the unitarity bound, 1 /2. The reason for this is not clear. But the lesson 
is that, at least in this simple model, to have a higher temperature superconductor, 
one should lower the dimension of the condensate. 



3.2 Conductivity 

We want to compute the optical conductivity, i.e. the conductivity as a function 
of frequency. By symmetry, it suffices to consider just the conductivity in the x 
direction. According to the gauge/gravity dictionary, this is obtained by solving for 
fluctuations in the Maxwell field in the bulk. Maxwell's equation with zero spatial 
momentum and time dependence e~ im gives the following equation for A x (r): 



We want to solve this with ingoing wave boundary conditions at the horizon, since 
this corresponds to causal propagation on the boundary, i.e., yields the retarded 
Green's function ll55ll . Asymptotically, 
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A,=Ai 0) + ^ + ... (17) 
r 

The gauge/gravity dictionary says the limit of the electric field in the bulk is the 
electric field on the boundary: E x = — A x °\ and the expectation value of the induced 
current is the first subleading term: J x = A X ■ From Ohm's law we get: 

, \ Jx Jx iJx iA* ^ , i n\ 

a{ - a) -E.—®~~^~-~^- (18) 

The real part of the conductivity is given in Fig. 3 for the case A = 2. Above 
the critical temperature, the conductivity is constant [33 1. As you start to lower the 
temperature below T c a gap opens up at low frequency. When these curves were first 
obtained, it was thought that Re[ff] was given by e~ A l T at small CO, which is what 
one expects from a BCS type description with an energy gap A . This would imply 
that at T = 0, Re [a] should vanish inside the gap. We will see that this is not the 
case. 
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Y 

Fig. 3 The formation of a gap in the real part of the conductivity as the temperature is lowered 
below the critical temperature. The curves describe successively lower temperatures. There is also 
a delta function at a> = 0. Figure is for the dimension two condensate and is taken from \ 29\. 




There is also a delta function at co = for all T < T c . This cannot be seen from a 
numerical solution of the real part, but it can be seen by looking for a pole in Im[<r]. 
A simple argument for this comes from the Drude model of a conductor. Suppose 
we have charge carriers with mass m, charge e, and number density n in a normal 
conductor. They satisfy 

dv v 
m — =eE — m— (19) 
dt T 

where the last term is a damping term and T is the relaxation time due to scattering. 
The current is J = env, so if E(t ) = Ee~' m , the conductivity is 
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kx 

°W - y— (20) 

where k — ne 2 /m. So 

kx kQyz 2 

Re[ff] = Im[(7] = (21) 

L J 1+co 2 t 2 L J 1+0) 2 T 2 

For superconductors, T — > °°, so Re[ff] oc 5(o) and Im[cr] « 1 /<a. 

A more general derivation comes from the Kramers -Kronig relations. These re- 
late the real and imaginary parts of any causal quantity, such as the conductivity, 
when expressed in frequency space. One of the relations is 

T r , 1 ^ r Re[a(o)')]do)' 

Imcr(fi)) = — & — — . (22) 

% J-oo co — CO 

From this formula we can see that the real part of the conductivity contains a delta 
function, if and only if the imaginary part has a pole. One finds that there is indeed 
a pole in Im[a] at co = for all T <T C . 

Fig. 4 shows the low temperature limit of the optical conductivity. The solid line 
denotes the real part and the dashed line denotes the imaginary part. The pole at 
CO = is clearly visible. 



1.5 r- 
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cr 
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Fig. 4 The low temperature limit of the optical conductivity for the dimension two condensate. 
The solid line denotes the real part and the dashed line denotes the imaginary part. Figure is taken 
from 03. 



A finite conductivity implies dissipation. In the bulk, this is reflected by the ingo- 
ing wave boundary conditions at the horizon. A standard normalizable perturbation 
of the Maxwell field would decay and get swallowed by the black hole. We main- 
tain a constant amplitude and purely harmonic time dependence by driving the mode 
with an applied electric field at the boundary. Dissipation normally causes a system 
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to heat up, and indeed if we had included the backreaction of the Maxwell pertur- 
bation on the metric, the flow of energy into the horizon would cause the black hole 
to grow and increase its temperature. 

When you approach the BF bound something interesting happens. As shown in 
Fig. 5, a new spike appears inside the gap. This looks like a bound state of quasi- 
particles with the binding energy given by distance between the pole and the edge 
of the gap. However, one must keep in mind that quasiparticles are a weak coupling 
concept and may not be well defined at strong coupling. We will have more to say 
about this spike in the next section. 
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Fig. 5 The low temperature limit of the conductivity for the dimension 3/2 condensate. Note the 
extra spike that appears inside the gap. Figure is taken from 1371 . 



Although we have focussed on the case of 2 + 1 dimensional superconductors 
with a 3 + 1 dimensional bulk dual, everything we have done is easily generalized 
to higher dimensions |37 1. The results are qualitatively the same. 

A robust feature that holds in both 2 + 1 and 3 + 1 superconductors and all X > 
Xbf is that 




with deviations of less than 10%JjThis is more than twice the weakly coupled BCS 
value of 3.5. Remarkably, measurements of this ratio in the high T c superconductors 
give roughly this value lfT8ll ! 



A few caveats should be noted: If one adds scalar interactions in the bulk by introducing a more 
general potential V(f ), the gap in the low temperature optical conductivity can become much less 
pronounced, so that this ratio becomes ill defined 1 23 1 . We will see in the next section that it also 
becomes ill defined at small q. Even when it is well defined, it is modified by higher order curvature 
corrections in bulk 1 1 9 1 . 
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4 Full solution with backreaction 

We now discuss the full solution to action (|2]i including the backreaction on the 
spacetime geometry. We start with the following ansatz for the metric 

ds 2 = -g(r)e- x ^dt 2 + —-+ r 2 (dx 2 + dy 2 ) (24) 

gin 

and use the same ansatz for the matter fields as before: 

A = ${r)dt, f=y/(7) (25) 
The equations of motion ar^J 

V' + ( 8 --^ + -)w' + q -^W--W = (26) 
\g 2 rj g l g 

*" + (y + ;)*'- 2 T^^° (27 > 

g 1 

--j)g+— 4 3r+^=0 (29) 

The first two are the scalar and Maxwell equations as before. The last two are the 
two independent components of Einstein's equation. (There are three nonzero com- 
ponents of Einstein's equation but only two are independent.) Note that the equa- 
tions for g and % are first order, and % is monotonically decreasing. These equations 



are invariant under a scaling symmetry analogous to ( 13 i: 



r^ar, {t,x,y) (t,x,y)/a, g^a 2 g, <j> -> a<$> . (30) 

When the horizon is at nonzero r, this can be used to set ro = 1 . These equations are 
also invariant under 

e x ^a 2 e x , t^at, Q-tQ/a (31) 

This symmetry can be used to set % — at the boundary at infinity, so the metric 
takes the standard AdS form asymptotically. 

The qualitative behavior of the solutions is similar to the probe limit, i.e., the 
probe limit indeed captures most of the physics. However there are a few important 
differences |30|. First, the apparent divergence in the dimension one condensate at 
low T is gone. Fig. [6] shows a plot of qO\ as a function of temperature. We have 
multiplied the condensate by the charge q since this is the quantity which is repre- 



9 We will set L = 1 for the rest of our discussion. 
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sented by the probe limit at large q. For q of order one, the curves are similar to the 
other condensates in the probe limit. But for large q, they grow at low temperature. 
One can show that for all q, the condensate still approaches a finite limit at T = 0. 
It is not clear why the dimension one condensate behaves differently in the probe 
limit from other masses. 



- 




0.0 0.2 0.4 0.6 0.8 1.0 
T 



Fig. 6 From bottom to top, the various curves correspond to q = 1, 3, 6, and 12. Figure is taken 
from |30). 



The second difference is more surprising. As explained earlier, the origin of the 
instability responsible for the scalar hair is the coupling of the charged scalar to 
the charge of the black hole. It was therefore expected that as q — > the instability 
would turn off. This is not what happens. A nearly extremal Reissner-Nordstrom 
AdS black hole remains unstable to forming neutral scalar hair, provided that m 2 
is close to the Breitenlohner-Freedman (BF) bound. This means that there is a new 
source of instability which can be understood as follows: An extremal Reissner- 
Nordstrom AdS black hole has a near horizon geometry AdS2 x R 2 . The BF bound 
for AdSd+i is m\ F = —d 2 /4. So scalars which are slightly above the BF bound for 
AdS4, can be below the bound for AdS2- This instability to forming neutral scalar 
hair is not associated with superconductivity (or superfluidity) since it doesn't break 
a t/(l) symmetry. At most it breaks a Z2 symmetry corresponding to \j/ — > — y/. Its 
interpretation in the dual field theory is not clear. 

One can make the following general argument for when an extremal Reissner- 
Nordstrom AdS black hole will be unstable to forming scalar hair lfl2l . Consider a 
scalar field with mass m and charge q in the near horizon geometry of an extremal 
Reissner-Nordstrom AdS black hole. Its field equation reduces to a wave equation 
in AdS2 with effective mass rn^ff = (>n 2 — 2q 2 )/6. The —2q 2 is the usual coupling 
to the electric charge and the factor of six comes from the difference between the 
radius of curvature in AdS<x and AdSx- The instability to form scalar hair at low 
temperature is then just the instability of scalar fields below the BF bound for AdS2- 
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m BF = ~ 1/4- Thus the condition for instability is 



2q 2 < -3/2 



(32) 



Of course, the mass must be above the four-dimensional BF bound, m 2 > —9/4, so 



that the asymptotic A1JS4 geometry remains stable. Eq. ( 32 1 is a sufficient condition 
to guarantee instability, but it is not necessary. 

To compute the conductivity, we again perturb the Maxwell field in the bulk. 
Assuming zero momentum and harmonic time dependence, the perturbed Maxwell 
field A x now couples to the perturbed metric component g tx . Physically, this is what 
one expects from the standpoint of the dual field theory. We are applying an elec- 
tric field and inducing a current. The current carries momentum, so T tx should be 
nonzero. This requires that the metric perturbation g tx must also be nonzercp*] How- 
ever, this also means that one has to solve for the thermal conductivity at the same 
time as the electrical conductivity. Fortunately, for homogeneous perturbations with 
harmonic time dependence, one can solve for g tx in terms of A x and get 





A' x + 


'{CD 2 <t>' 2 \ x 2q 2 Y 2 ' 


.8 2_ 




\g 2 8 ) 8 



0. 



(33) 



This is similar to ( 16 1 except for the form of the metric and the extra term <j> /g 



coming from the metric perturbation g tx . The conductivity is still 



(7(0)) 



I) 



i A x 



(34) 



The results for the optical conductivity are qualitatively similar to the probe limit, 
and shown in Fig. 7. There are a few differences. One is that the conductivity in the 
normal phase is no longer constant. Another difference is that in the low temper- 
ature limit, the gap in Rec(a)) for small 0) becomes less pronounced at small q. 
The robust feature (O g /T c w 8 seen in the probe limit continues to hold for q > 3, 
but is less robust for q < 3 mainly because the gap becomes less well defined. Per- 
haps the most important difference is that Red now has a delta function contribution 
at (0 — even in the normal phase when T > T c . This infinite DC conductivity is 
not superconductivity, but just a consequence of translational invariance. A transla- 
tionally invariant, charged system does not have a finite DC conductivity because 
application of an electric field will cause uniform acceleration. One does not see 
this in the probe limit, since we fixed the gravitational background which implicitly 
breaks translation invariance. One indication of this is that we get an electric current 
without momentum flow as we mentioned earlier. Although there is a delta function 
at co = for all temperature, its coefficient grows as the temperature is lowered be- 
low T c . This indicates the presence of a new contribution coming from the onset of 
superconductivity. 



' I thank Hong Liu for this comment. 
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CO 

T 

Fig. 7 Conductivity for the dimension two condensate with q = 3. The dashed line is the real part 
of the conductivity at T = T c and the solid lines are the same conductivities at successively lower 
temperature T/T c = 0.651, 0.304. There is a delta function at the origin in all cases. Figure is taken 
from |30). 



4.1 Reformulation of the conductivity 



Equation ( 33 1 can be simplified by introducing a new radial variable 



e X/2 

dz= dr (35) 



At large r, dz = dr/r 2 , and we can choose the additive constant so that z = — 1 jr + 
0(1 /r 2 ). Since g vanishes at least linearly at a horizon and % is monotonically 



decreasing, the horizon corresponds to z = — °°. In terms of z, (33 I takes the form of 
a standard Schrodinger equation with energy (O 2 : 

-A XJZ + V(z)A x = eo 2 A x (36) 

where 

V{z)=g[<t> 2 r + 2q 2 y 2 e-X] (37) 

From the known asymptotic behavior of the solution near infinity, one can show 
that V(0) = if the dimension of the condensate, X, is greater than one, V(0) is a 
nonzero constant if A = 1, and V(z) diverges asz— >0ifl/2<A < 1. One can also 
show that the potential always vanishes at the horizon l38ll . 



We want to solve ( 36 1 with boundary conditions at z = — °° corresponding to 
waves propagating to the left. This corresponds to ingoing boundary conditions at 
the horizon which is needed to extract causal results. The easiest way to do this is 
to first extend the definition of the potential to all z by setting V = for z > 0. An 
incoming wave from the right will be partly transmitted and partly reflected by the 
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potential barrier. Since the transmitted wave is purely ingoing at the horizon, this 
satisfies our desired boundary conditions. Writing the solution for z > as A x = 

e~ i(0z + Me imz , we clearly have A x (0) = 1 +M and A x . z {0) = -i(o(l -M). In terms 



of z , A x 1] = -A X)Z (0), so from <|34i 



a (a) = - — -= (38) 

The conductivity is directly related to the reflection coefficient, with the frequency 
simply giving the incident energy [38J! The qualitative behavior of (7(a)) is now 
clear. Let us first assume that the dimension of the condensate is X > 1 so that 
V is bounded. At frequencies below the height of the barrier, the probability of 
transmission will be small, Iffl will be close to one, and o((0) will be small. At 
frequencies above the height of the barrier, 3fc will be very small and o(co) ~ 1 (the 
normal state value). Clearly the size of the gap in o((o) is set by the height of the 
barrier: co g ~ \fV mux . The case 1 /2 < A < 1 is qualitatively similar. Even though 
the potential is not bounded, y/V is integrable, so there is still tunneling through the 
barrier. 

As one lowers the temperature, the condensate yf increases and the potential be- 
comes both higher and wider (see Fig. |HJ. This causes an increasing exponential 
suppression which was erroneously interpreted as evidence for Re ff~ e~ A l T in the 
original papers on holographic superconductors. However, as T — > 0, the potential 
approaches a finite, limiting form which still vanishes at z = — °°. This means that 
there will always be nonzero tunneling probability and hence a nonzero conductivity 
even at small frequency. In other words, there is no hard gap in the optical conduc- 
tivity at T = 0. This conclusion continues to hold if the m 2 |f| 2 term in our bulk 
lagrangian is replaced by a general potential V(f ). The fact that there is not a hard 
gap is consistent with calculations of the specific heat which show that it vanishes 
as a power law at low temperature and not exponentially. 

As discussed in section 3.2, in order for Rec to have a delta function at CO = 
representing the infinite DC conductivity, one needs Im a to have a pole at CO = 0. It 
is easy to see that this is indeed the case, for any positive Schrodinger potential V(z) 



that vanishes at z = — °°. Imagine solving ( 36 1 with (0 = 0, and A x = 1 at z = — °°. 
(This represents the normalizable solution and is the zero frequency limit of the 
ingoing wave boundary condition.) Since A x > 0, the solution will be monoton- 

ically increasing. At z = 0, A* = A X (Q) and A* = — A x>z (0). These are both real 
and nonzero. From p4| ) it then follows that Imff has a pole at CO = 0. Even in the 
normal phase when yf = 0, the potential is nonzero due to the contribution from the 
electric field. This is precisely the term which arises due to the backreaction on the 
metric and is absent in the probe approximation. This is another way to see why the 
probe limit does not have a delta function at co = in the normal phase while the 
full backreacted solution does. Note that the Schrodinger potential remains finite in 
the probe limit since qyi is held fixed. 

This approach also explains the spike in the conductivity that was seen in Fig. 5. 
At low frequency, the incoming wave from the right is almost entirely reflected. If 
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Fig. 8 Schrodinger potential for X = 2, q = 10. The top curve as T = and the bottom curve has 
T = T c . Figure is taken from 1 38 1. 



the potential is high enough, one can raise the frequency so that about one wave- 
length fits between the potential and z = 0. In this case, the reflected wave can 
interfere destructively with the incident wave and cause its amplitude at z = to be 
exponentially small. This produces a spike in the conductivity. If one can raise the 
frequency so that two wavelengths fit between the potential and z = one gets a 
second spike, etc. More precisely, using standard WKB matching formula, spikes 
will occur when there exists ft) satisfying 



ft) z 



-V(z)dz- 



K 

4 



n% 



(39) 



for some integer «, where V(— zo) = ft) 2 . It is now clear that the spike which we 
saw in the probe approximation with m 2 saturating the BF bound, also appears in 
the full backreacted solution, and for some m 2 slightly above the BF bound. When 
the spikes were first seen, it was speculated that they corresponded to vector normal 
modes of the hairy black hole. It is now clear that they are not true normal modes 
even at T = 0, since A x does not actually vanish at infinity. The actual modes all 
have complex frequency and correspond to familiar quasinormal modes. In other 
words, there are no bound states in this potential (with boundary condition A x = 
at z = 0) since the potential vanishes at z = — °°. 



5 Zero temperature limit 

The extremal Reissner-Nordstrom AdS black hole has large entropy at T — 0. If this 
was dual to a condensed matter system, it would mean the ground state was highly 
degenerate. The extremal limit of the hairy black holes dual to the superconductor is 
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not like Reissner-Nordstrom. It has zero horizon area, consistent with a nondegen- 
erate ground state ||38l . One might wonder how this could be the limit of the T > 
black holes since we have seen that for all T ^ 0, one can scale r so that ro = 1. The 
point is simply that the dimensionless ratio ro/^^Oas7'^0. 

The extremal limit also has zero charge inside the horizon (except in the case 
q = when the scalar hair cannot carry the charge). This is an expected conse- 
quence of the horizon having zero area. If the black hole tried to keep a nonzero 
charge as the temperature was lowered and its horizon shrunk to zero, the electric 
field on the surface would grow arbitrarily large. Eventually this would pair create 
charged particles via the Schwinger mechanism and neutralize the black hole. This 
is a quantum argument involving pair creation, but it has a classical analog in terms 
of a superradiant instability for charged scalar fields. 

The near horizon behavior of the zero temperature solution depends on the mass 
and charge of the bulk scalar field. We will discuss two classes of solutions below. 
In both cases, one can solve for the solutions analytically near r — 0. These leading 
order solutions depend on a free parameter which can be adjusted so that the solution 
asymptotically satisfies the desired boundary condition. In the cases we study, the 
solution is typically not smooth at r — 0, but since they are the limit of smooth black 
holes, they are physically allowed. If one modifies the potential for f so that it has 
more than one extremum, then smooth, zero temperature solutions exist in which f 
rolls from one extremum to another [23]. 



5.1 m 2 = 



This corresponds to a marginal, dimension three operator developing a nonzero ex- 
pectation value in the dual superconductor. To determine the leading order behavior 
near r = 0, we try an ansatz 



(40) 

We have used the scaling symmetries ( |30] > and ( 3 1 1 to set the coefficients in (j) and 
g to one. Substituting this into the field equations and equating the dominant terms 
for small r (assuming a > —1), one finds: 



q¥o = 
a + 2 



a- + 5a + 6 
2 



1/2 



Xi = 



qe 



Xo 



4(a + l) 



cr + 5a + 6 



1/2 



(41) 



(42) 



2(2a 2 + 7a + 5) 

One can now numerically integrate this solution to large radius and adjust a 
so that the solution for y is normalizable. One finds that this is possible provided 
q 2 > 3/4. This is consistent with the condition for instability (32i. The value of a 
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depends weakly on q (see Fig. 9). In all cases, \a\ < .3. Fig. 10 shows the results 
for y/ and g for the zero temperature solution, and shows how the T > solutions 
approach it as T 0. 



0.3 




2 

q 

Fig. 9 Values of a for various charges. Note that it approaches a constant as <y — > oo. Figure taken 
from OS). 
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Fig. 10 Zero temperature, A = 3 and q = 1 solution (dashed line), compared to successively lower 
temperature hairy black holes (solid lines.) Note that g almost has a double zero at r/fl = 1 /2\/3 
which is the extremal horizon for Reissner-Nordstrom AdS. In the limit q — > \/3/2 the solution 
becomes extremal Reissner-Nordstrom AdS with all hair behind the horizon. Figure is taken from 

EH). 



Near r = 0, % approaches a constant and g~r 2 . Thus the metric approaches AdS4 
with the same value of the cosmological constant as the asymptotic region. The 
extremal horizon is just the Poincare horizon of AdS a ■ The scalar field approaches a 
constant and the Maxwell field vanishes. In terms of the dual field theory, this means 
that the full conformal symmetry is restored in the infrared. 
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These solutions are not singular at r = since all curvature invariants remain 
finite. They can be viewed as static, charged domain walls. Even though they are not 
singular, when a ^ the solutions are not C°° across r — 0. Some derivatives of the 
curvature will blow up. However, there is a special value of the charge, q = 1.018, 
where a = 0. This solution is completely smooth across the horizon. It is not clear 
what the significance of this value is for the dual field theory. 



5.2 q 2 > \m 2 \/6 (m 2 < 0) 



We try the following ansatz near r = 0: 

V/ = A(-logr) 1/2 , g = g r 2 (- log r), (j) = (j) Q r P (-logr) 



1/2 



(43) 



where we have used the radial scaling symmetry ( 30 1 to set an arbitrary length scale 



in the logarithm to one. The behavior of % is determined by ( 28 1 and whether we 
expect ry/' 2 or rq 2 ^ 2 yf 2 e x / g 2 to dominate. We assume ry/' 2 dominates, so 



e x = K(~logr) 



A 2 / 4 



(44) 



with K a constant of integration. The second term in ( 28 1 is indeed negligible pro- 
vided j3 > 1. Equating the dominant terms in the equations of motion leads to: 



3 



1 ±- [ 1 

2 2 



48^ 



2\ 1/2 



(45) 



It is clear that for appropriate metric signature we need positive go which tells 
us this ansatz is only appropriate for negative m 2 . Since we require J3 > 1, only the 



plus sign in (45 i is allowed and we require 

q 2 > -m 2 /6 



(46) 



With these restrictions, our near horizon solution is 

y/ = 2(-logr) 1/2 , g = (2m 2 /3)r 2 logr, e x = -Klogr (47) 

</) = ^ /3 (-logr) 1 / 2 , (48) 

The one remaining free parameter is 0o- This parameter can be tuned so that the 
asymptotic boundary condition on y/ is satisfied. 

The horizon at r = has a mild singularity. The scalar field diverges logarithmi- 
cally and the metric takes the form (after rescaling t) 



ds 2 = r 2 {—dt 2 + dxjdx') - 



dr 2 
g Q r 2 (-logr) 



(49) 
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Notice that Poincare invariance is restored near the horizon, but not the full con- 
formal invariance. (Some early indications of emergent Poincare symmetry were 
found in ll23ll ). Introducing a new radial coordinate, r = — 2(— logr) 1//2 /gV 2 the 
metric becomes 

ds 2 = e- 8of2/2 (-dt 2 + dxidx') + df 2 (50) 
near the horizon which is now at r = — °°. 



6 Adding magnetic fields 

As we mentioned in the introduction, one of the characteristic properties of super- 
conductors is that they expel magnetic fields. However, a superconductor can expel 
a magnetic field only up to a point. A sufficiently strong field will destroy the super- 
conductivity. One defines the thermodynamic critical field B c by setting the energy 
it takes to expel the magnetic field equal to the difference in free energy between the 
normal and superconducting states: 

&^=F„(T)-F S (T) (51) 

67Z 

where V is the volume. Superconductors are divided into two classes depending 
on how they make the transition from a superconducting to a normal state as the 
magnetic field is increased. In type I superconductors, there is a first order phase 
transition at B = B c , above which magnetic field lines penetrate uniformly and the 
material no longer superconducts. In type II superconductors, there is a more gradual 
second order phase transition. The magnetic field starts to penetrate the supercon- 
ductor in the form of vortices with quantized flux when B = B c \ < B c . The vortices 
become more dense as the magnetic field is increased, and at an upper critical field 
strength, B = B C 2 > B c , the material ceases to superconduct. 

We now show that our 2+1 dimensional holographic superconductors are type II. 
(Interestingly enough, the high T c superconductors are also type II.) We will view 
this superconductor as a thin film superconductor in 3+1 dimensions with a perpen- 
dicular magnetic field. Suppose the 2+1 dimensional sample is a disk of radius R. 
In order for the disk to expel the magnetic field, it must produce a current circu- 
lating around the perimeter. Solving Maxwell's equations, this current will expel a 
field not only in the area kR 2 of the disk but in a larger volume of size V3 ~ R 3 . 
Since the difference in free energies between the normal and superconducting state 
scales like R 2 , in the large R (thermodynamic) limit, the superconductor does not 
have enough free energy available to expel a magnetic field. Thus magnetic fields of 
any non-vanishing strength will penetrate a thin superconducting film and B c = 0. 

To show that the holographic superconductor is type II, we start in the normal 
phase with a large applied magnetic field, and slowly lower B. We will see that parts 
of the system start to become superconducting at a value B C 2 > lfTl [30ll . To model 
the normal phase with an applied magnetic field, we use a black hole in AdS with 
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both electric and magnetic charges (but no scalar hair). The magnetic field in the 
bulk corresponds to applying a background magnetic field on the boundary. The 



metric takes the form ( 24 1 with % = and 



g(r) = r 2 -^ (4r* + P 2 +B 2 ) + ^r 2 (p 2 +B 2 ) . (52) 
The black hole temperature is 

T = — . (53) 



The vector potential is 



A=p ( ) at B.xdy. (54) 



Since we are interested in the onset of superconductivity which is a second order 
phase transition, the condensate will be small and we can ignore the backreaction of 
y/ on the metric. (Note that this is different from the probe approximation since the 
backreaction of the Maxwell field is included here.) Since the gauge potential now 
depends on x, we separate variables as follows 

Y(r,x,y)=R(r)X(x)eto (55) 

The equation for X{x) turns out to be a simple harmonic oscillator centered at k/qB 
with width (qB)~ l l 2 . The eigenvalues are 2qB(n + i) but only the lowest mode is 
expected to be stable. Note that the momentum k shifts the origin of the harmonic 
oscillator, but does not change its energy. This is just the degeneracy of Landau 
levels. The radial equation becomes: 



(r 2 8 R'y 



1 2 P 2 ( r A , „ D .„2„2 



- 1 + qB-m z r / 



R = (56) 



One can now choose B so that there is a normalizable solution. This defines the up- 
per critical field B c 2- One finds that B C 2 is a decreasing function of temperature and 
vanishes at T = T c . This is what one expects: At lower temperatures, a supercon- 
ductor can support a larger magnetic field. Due to the degeneracy in k, the actual 
condensate involves a superposition of these modes for all k and forms a lattice of 
vortices. Maeda et al. Il46l have recently studied the free energy as a function of 
the two parameters which govern the shape of the vortex array and shown that the 
minimum of the free energy at long wavelength corresponds to a triangular array 
which is what is predicted by the Landau-Ginzburg model. 
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As mentioned in section 2, the boundary theory does not have dynamical Maxwell 
fields, so one cannot actually see the Meisner effect. However, one can show that in 
the superconducting state, static magnetic fields induce currents whose backreaction 
would cancel the applied magnetic field B31 [30l . In fact, one can show that these 
currents are proportional to A, and thus reproduce London's equation. 

To make life easier, in this section we shall work in the probe limit (q — > °°) in 
which the metric is kept fixed to be simply the Schwarzschild AdS black hole ([3]). 
Assume that we have solved for A, and y/ in this background as described in section 
(3.1). We then wish to add a static perturbation A x that has nonzero momentum. 
We will take the momentum to lie in a direction orthogonal to A^, so we can con- 
sistently perturb the gauge field without sourcing any other fields. This perturbation 
corresponds to adding a static magnetic field in the bulk, and the asymptotic value of 
this magnetic field corresponds to the magnetic field added to the boundary theory. 
As before, the first subleading term in A x (r) corresponds to the induced current. 

Rather than doing this calculation explicitly, we can obtain the answer by relat- 
ing it to the calculation in section 3.2. Let A x have both momentum and frequency 
dependence of the form A v ~ e -uot+iky ^ e ra( jj a i equation for A x (r) reduces to 

CO 2 k 2 \ . 



(fK) + ( y - ^ ) A* = VA* , (57) 

where ' denotes differentiation with respect to r. We studied this equation with k = 
in section 3.2, and found a pole in Im(7[o] at CO — 0. Let the residue of this pole be 



n s , the superfluid density. From dl8j> this means that A* = —n s A x °K But J x = A X 



and A x °^ is just the vector potential in the boundary theory, so 

J x {co,k = 0) = -n s A x (d),k = 0) (58) 
It is clear from the form of eq. p7| ) that the limits co — > and k — > must commute. 



To compute n s , we simply set both co and k to zero and solve (57 1. Thus for small k, 
we directly obtain the London equation 

Ji(a>,k) = -n s Ai(a>,k) (59) 

where we have used the isotropy of space. This equation is clearly not gauge invari- 
ant. It is expected to hold in a gauge where V;A' = 0. 



6.2 Correlation length 



The gauge/gravity dictionary says that the retarded Green's function (for J x ) in the 
dual field theory is 
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A V° Kk [ (60) 
AP((0,k) 

We can define a correlation length S, by expanding 

ReG R (0,k) = -n s (l+% 2 k 2 + ---) (61) 



By solving ( 57 1 numerically, including the k dependence, one can compute this 
correlation length. One finds that it diverges near the critical temperature like 
% T c w 0. 1 ( 1 - T I T c ) - 1 1 2 El . A similar diverg ence is found in a correlation length 
obtained from fluctuations in the condensate [45|. 



6.3 Vortices 

As we discussed above, at the onset of superconductivity when B = B c i, there is a 
lattice of vortices. However at this point, the condensate is small everywhere and the 
equations could be linearized. It is of interest to find the bulk solution describing a 
vortex lattice away from B = B c o. The first step is to find the single vortex solution. 
This corresponds to a something like a cosmic string stretching from the black hole 
horizon to infinity. However there is a crucial difference between a cosmic string 
and the dual of a vortex. A standard cosmic string has a fixed proper radius. This 



means that its size in the x, y coordinates used in ( 24 1 goes to zero at large r. In terms 
of the boundary theory, this is a point-like object. We want a solution in which the 
cosmic string is a fixed size in the x,y coordinates, so it has a finite radius in the 
superconductor. This looks more like a funnel in the bulk with a proper radius that 
grows as one increases r. 

Such a solution has recently been found in the probe limit [49 |{^]Let us write 
the background metric using polar coordinates for the flat transverse space: 

dr 2 

ds 2 = -f(r)dt 2 + ■— r + r 2 (dr 2 + r 2 d(p 2 ) . (62) 

Assume an ansatz 

W=lif(r,f)e in 1', A r =A r (r,r), A^^A^r) (63) 

Substituting this into the field equations, one obtains a set of nonlinear PDE's. The 
condensate is now a function of radius r. It vanishes at r = which represents the 
center of the vortex, and approaches a constant at large r. This shows that there is no 
superconductivity inside the vortex. Since the Maxwell field is not dynamical on the 
boundary, the magnetic field is not localized inside the vortex as one would expect. 
Instead, the authors of ll49l work in finite volume and assume a uniform magnetic 
field on the boundary with total flux equal to 2nn as expected for the vortex. 



1 1 See 1 2 1 for another approach. 
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7 Recent developments 

So far, we have discussed planar black holes in the bulk. Spherical black holes have 
the same instability and describe spherical superconductors. One can add rotation 
to a spherical black hole, and the effect of this rotation on the dual superconductor 
has recently been studied [56]. It was found that for T < T c (the transition temper- 
ature at zero rotation), there is a critical value of the rotation which destroys the 
superconductivity in analogy to the critical magnetic field. 

The bulk gravitational theory that we have studied was not derived from string 
theory. It was just the simplest model that could describe a dual superconductor, and 
we have seen that it works rather well. It predicts that a charged condensate forms 
at low temperature (in a second order phase transition), that the DC conductivity 
is infinite, and that there is a gap in the optical conductivity at low frequency - all 
basic properties of superconductors. However, to go beyond this and have a more 
detailed microscopic understanding of the dual superconductor, one needs to embed 
this model in string theory. A linearized instability of the type needed to describe 
a holographic superconductor was found in a string compactification in [12|. This 
has now been extended to a full description by two different groups. Gauntlett et al. 
ifTBI realized the m 2 — —2, q = 2 model in M theory. In other words, they found a 
consistent truncation of eleven dimensional supergravity in which the four dimen- 
sional fields were just a metric, Maxwell field and charged scalar with this mass 
and charg^] Gubser et al. Il26l realized the same model in one higher dimension 
(a five dimensional bulk which is dual to a 3 + 1 dimensional superconductor) with 
m 2 = —3, q = 2 in type IIB string theory. Both groups used Sasaki-Einstein com- 
pactifications with U(l) symmetry, where the charged scalar is related to the size of 
the U(l) fibration. 

There are two main differences between the models obtained from string theory 
and our simple model. The first is that the kinetic term for the charged scalar in- 
volves nonminimal coupling. (This was investigated earlier by Franco et al. Q3).) 
The second difference is that the scalar potentials V(y) coming from supergravity 
are more complicated than the simple mass terms we have considered so far. This 
does not affect the onset of superconductivity: since the scalar field is small near T c , 
the critical temperature is just determined by the mass term in the potential. How- 
ever, the low temperature limit is different. The potentials coming from supergravity 
have more than one extremum. In this case, the zero temperature limit of the hairy 
black hole is a smooth domain wall in which the scalar rolls from one extremum 
in the UV to another extremum in the IR ll24ll . The IR limit is another copy of AdS 
with, in general, a different cosmological constant. These examples all have emer- 
gent conformal symmetry in the IR. 

The optical conductivity in these theories never show a pronounced gap at low 
temperature. Typically, Rea(o)) = k(0 S for small co with a coefficient k of order 
one. From our Schrodinger equation reformulation this is a bit surprising. Why isn't 



In general, there is an additional neutral scalar, but for purely electrically charged black holes, 
this field can be set to zero. 
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there an exponential suppression due to tunneling through the barrier? The answer is 
simply that the effective potential in the Schrodinger equation never gets very large 
in these solutions since all fields are order one. 

Even though embedding in string theory makes it possible to study the dual field 
theory in more detail, this is still difficult in 2 + 1 dimensions. It is somewhat easier 
in 3 + 1 dimensions since the dual theories are related to ordinary gauge theory. One 
particularly interesting aspect of the 3 + 1 superconductor is that the condensate 
includes a term bilinear in the fermions, like a Cooper pair |26 1. 

It was shown in |22| that one can replace the Maxwell field and charged scalar 
in the bulk with an SI/ (2) gauge field and still have a dual description of a super- 
conductor. The main difference is the following: The symmetry of a superconductor 
refers to the energy gap above the Fermi surface. The holographic superconduc- 
tors we have discussed here all have S-wave symmetry since they do not prefer a 
direction. With an SU (2) gauge field in the bulk, one obtains a P-wave supercon- 
ductoip^| This leads to another way to realize holographic superconductors (in the 
probe limit) in string theory: one can use the SU (2) gauge field realized on a pair of 
branes flUED. 

This is a rapidly evolving field and there have been several developments since 
the Milos summer school in September 2009, including: 

• A key property of condensed matter systems is, of course, the atomic lattice. This 
gives rise to the phonons in the BCS theory which cause the electrons to pair. Our 
simple model has translational symmetry and no sign of a lattice. However, some 
interesting ideas on how to include a lattice have been discussed by Kachru et al. 

HO- 

• Some superconductors exhibit striped phases associated with charged density 
waves. It was shown in 1 50] that the Chern-Simons term in five dimensional 
supergravity can lead to an instability at nonzero momentum in which a spatially 
modulated condensate forms. 

• The response of fermions to the superconducting phase has been studied lfTTI IT4l 
[231 . It was found that with a suitable coupling between the bulk fermions and 
scalar, there are stable quasiparticles with a gap 1141 . 



8 Conclusions and open problems 

For ninety years after its discovery in 1915, general relativity was viewed as a the- 
ory of gravity. It has proved very successful in describing a wide range of gravita- 
tional phenomena from the bending of light to gravitational waves and black holes. 



One cannot compute the energy gap directly. Indeed, as we have seen, there probably is not a 
strict gap. The interpretation as a P-wave superconductor comes from the fact that there is a vector 
order parameter and the conductivity is strongly anisotropic in a manor consistent with P-wave 
nodes on the Fermi surface. 
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But in the last few years we have seen that this same theory can describe other 
areas of physics as well, including superconductivity. This is all due to the magic 
of gauge/gravity duality. Although the full power of gauge/gravity duality relates a 
quantum theory of gravity (indeed string theory) in the bulk to a nongravitational 
theory on the boundary, we have worked in a large N limit in which the bulk theory 
is just classical general relativity. This large N limit also explains how we can have 
spontaneous symmetry breaking in a 2 + 1 dimensional field theory, in apparent con- 
tradiction to the Coleman-Mermin-Wagner theorem. The large N limit suppresses 
fluctuations in the fields^ 

It is natural to ask how surprised one should be that general relativity can repro- 
duce the basic properties of superconductors. After all, Weinberg [60 1 has shown 
that much of the phenomenology of superconductivity follows just from the spon- 
taneous breaking of the 1/(1) symmetry. Once we have found the instability that 
leads to charged scalar hair, doesn't everything else follow? There are indications 
that something deeper is going on. For example, order one dimensionless ratios can 
be computed and compared with experiment. In particular, the ratio CO g /T c ss 8 dis- 
cussed in section 3.2, is close to the observed value. This does not follow from 
symmetry arguments alone. 

Since our bulk dual of a superconductor just involves gravity interacting with a 
Maxwell field and a charged scalar, there is a superficial similarity to a Landau- 
Ginzburg description. However, it is important to keep in mind two key differences. 
First, the low temperature instability must be put in by hand in the Landau-Ginzburg 
model, whereas it arises naturally in our gravitational description. Indeed, we have 
seen that there are two physically distinct instabilities which can trigger the phase 
transition. Second, the Landau-Ginzburg model is only valid near the transition tem- 
perature, since it involves a power series in the order parameter (p. To go beyond 
T « T c , one would need to specify an entire potential V((p). Initially, our bulk the- 
ory also had the freedom to add an arbitrary potential V(f ). We chose just a mass 
term for simplicity. However, once one embeds the bulk theory into string theory, 
the potential is fixed and is no longer arbitrary. 



8.1 Open problems 

We close with a list of open problems^ They are roughly ordered in difficulty with 
the easier problems listed first. (Of course, this is my subjective impression. With 
the right approach, an apparently difficult problem may become easy!) 

1. In the probe limit below the critical temperature, there is an infinite discrete set 
of solutions for y/ which are all regular on the horizon and satisfy the required 
boundary condition at infinity. (When backreaction is included, there is only a 



This is not an issue for the 3 + 1 dimensional superconductor, which can be holographically 
described by the bulk gravitational theory |2j in one higher dimension. 
15 I thank Sean Hartnoll for suggesting some of these problems. 
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finite number of such solutions.) They can be labeled by the number of nodes 
(zeros) they contain. It is widely believed that only the lowest solution with no 
nodes is stable. Is this true in the theory we have discussed where the scalar 
potential V (y/) just contains a mass term? Is it true if the theory contains scalar 
self interactions, like the potentials derived from string theory? 

2. We have seen that an extreme Reissner-Nordstrom AdS black hole is unstable to 
forming neutral scalar hair if the scalar has mass close to the BF bound ([9]). Under 
evolution, as the black hole is developing scalar hair, the total charge and mass 
at infinity are conserved. Since the scalar field can carry energy, but not charge, 
the final black hole (which need not be extremal) has the same charge, but less 
mass than the initial extremal limit. How is this consistent with the original black 
hole being extremal? What is the extremal limit in this Einstein-Maxwell-scalar 
theory? A related question is: what is the interpretation of this neutral scalar hair 
in the dual field theory? 

3. We studied magnetic fields in 2+1 dimensional superconductors in section 6. 
What happens if you add a magnetic field to 3+1 superconductors? One immedi- 
ate difference is that B c is no longer zero. Are the holographic superconductors 
still type II? The new magnetic brane solutions in 1 13] may be useful. 

4. As discussed in section 6.2, a vortex solution has recently been found in the probe 
approximation. Find the vortex solution with backreaction. More generally, find 
the solution describing a lattice of vortices in the bulk. This would correspond 
to a magnetically charged black hole in which the magnetic flux is confined to 
"cosmic strings" stretching from the horizon to infinity. To solve this problem, 
one must solve nonlinear partial differential equations. 

5. Recall that the symmetry of a superconductor refers to the energy gap on the 
Fermi surface. The holographic superconductors we have discussed here all have 
S-wave symmetry since they do not prefer a direction. P-wave holographic super- 
conductors have been found [22 54]. The high T c cuprates are known to have D- 
wave symmetry. Can one find a holographic superconductor with D-wave sym- 
metry? It is natural to try to condense a charged spin two field. Interestingly 
enough, the newly discovered iron pnictides appear to have S-wave symmetry, 
but it is complicated by multiple Fermi surfaces 111 Oil . 

6. It was discovered recently that there are materials in which the (DC) conductivity 
goes to zero for T < T c |58 1. This is the opposite of a superconductor and is called 
a "superinsulator". Can one find a holographic description of a superinsulator? 
It is tempting to use the fact that under electromagnetic duality in AdSn, the 
conductivity transforms as a — > — 1 /cr |28l . However one must start with a bulk 
action which is invariant under electromagnetic duality. The action (|2|i is not, 
since it explicitly involves A^. 

7. As we discussed in section 7, it has recently been shown that the holographic 
superconductors can be derived from string theory, so in some cases the dual mi- 
croscopic theory is known. Can one understand the pairing mechanism in these 
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cases (if indeed "pairing mechanism" is the right concept at strong coupling)? 
This remains one of the main open questions in understanding the high T c mate- 
rials. 

8. As mentioned above, the large N limit is responsible for allowing spontaneous 
symmetry breaking in our 2+1 dimensional field theory. Can one show that away 
from this limit, massless fluctuations lead to infrared diverges which destroy the 
long range order? 

9. In gauge/gravity duality, one takes a large N limit to justify using classical gen- 
eral relativity in the bulk. What is the analog of this large N limit in condensed 
matter systems? In other words, what types of materials are likely to have a 
(tractable) dual gravitational gravitational description? (See [52] for a discussion 
of some of the issues.) 

10. The high temperature cuprate superconductors satisfy a simple scaling law re- 
lating the superfluid density, the normal state (DC) conductivity and the critical 
temperature [36|. Can this be given a dual gravitational interpretation? 
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